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Abstract 

A fluid dynamics video is linked to this article, which have been 
submitted to the Gallery of Fluid Motion as part of the 65th American 
Physical Society meeting of the Division of Fluid Dynamics, held in 
San Diego, Cahfornia, USA, over 17-20 November 2012. The video 
serves to visualize flows generated in a rectangular enclosure that are 
subjected to both mechanical and thermal forcing through a common 
horizontal boundary. This system exhibits features consistent with 
either horizontal convection or lid-driven cavity flows depending on 
the ratio between thermal and mechanical stirring, and three different 
cases are visualized in the linked videos. 

1 Introduction and video links 

The video contained in this submission depicts a combination of thermal 
forcing through a horizontal convection mechanism [31 [2] and mechanical 
forcing through a lid-driven mechanism [3] being applied to a quiescent 
body of fluid contained within a rectangular enclosure. In the video, both 
thermal and mechanical forcing are applied through the bottom boundary. 
Hence the depicted flow configuration combines both horizontal convection 
dynamics and lid-driven cavity dynamics, and may serve as a model for 
global overturning ocean current systems such as the Meridional Overturn- 



ing Circulation (MOC) [6j. Both a high-resolution version of the video (720i^ 
|19.0 MB)| and a [lower-resolution version (6.97 MB)| are available. 



2 Methodology 

The flow sequences are obtained numerically using an incompressible spectral- 
element solver combined with a Boussinesq treatment of buoyancy [5j . The 
results are also representative of configurations where forcing is imposed 
along the top boundary, with thermal forcing being inverted (because cooler 
fluid descends with the same dynamics as hotter fluid rises). 

The flows are visualized by plotting contours of the temperature field, 
which are proportional to buoyancy and inversely proportional to density 
in the Boussinesq approximation. Variations in the temperature field are 
emphasized by slightly elevating the visualization plane in the out-of-plane 
direction according to temperature, and lighting the surface elevation to 
produce a three-dimensional effect. 

The enclosure has a ratio of height to width of 0.16, and the width L 
serves as the reference length for definitions of both the Rayleigh number 
and Reynolds number. The Rayleigh number characterizes the strength of 
the thermal forcing applied as a linearly varying temperature distribution 
along the bottom boundary, and is defined as 

Ha = , 

where g is the gravitational acceleration, 6T is the temperature difference 
imposed along the forcing boundary, i' is the fluid kinematic viscosity, a 
is the fluid volumetric expansion coefficient, and ht is the ffuid thermal 
diffusivity. The Prandtl number of the fluid is defined as 

V 

Pr= — , 

and is fixed at Pr = 6.14 for all three cases in this video, which corresponds 
to water at standard laboratory conditions. 

The Reynolds number characterizes the strength of mechanical stirring 
applied to the fluid through a constant horizontal motion of the bottom 
boundary from left to right (in the direction of increasing heating, and there- 
fore increasing added buoyancy). It is deflned as 

V 

where V is the velocity of the moving boundary. 
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3 Video contents 



Three cases are depicted. These sequences, which each run for approxi- 
mately 40 seconds, show the transient startup phase of the flow under dif- 
ferent conditions as follows: 



1. 


Ra = 


10^°, Re = 


4072 


2. 


Ra = 


10^°, Re = 


5 700 


3. 


Ra = 


10^°, Re = 


7329 



The common high Rayleigh number ensures that in the absence of any 
mechanical stirring (i.e. Re = 0) the flow is in a convection-dominated 
regime. At this condition the long-timescale equilibrium state features con- 
tinuous time-dependent eruption of buoyant plumes of hot fluid at the heated 
right-hand end of the enclosure, driving a counter-clockwise overturning cir- 
culation in the enclosure that is completed by a diffusive return of fluid at 
the cooler left end of the enclosure. The three Reynolds numbers straddle a 
zone in which a gradual shift occurs from a flow with overturning circulation 
dominated by natural convection (Case 1), through a mixed regime (Case 2) 
to a flow primarily exhibiting features associated with forced convection and 
mechanical base-driven stirring (Case 3). 

An explanation for the different behaviors observed in each of these se- 
quences can be determined by considering the velocity boundary layers pro- 
duced at the forcing boundary by each of these mechanisms. In a convection- 
dominated flow (e.g. Case 1), the buoyancy destabilization invokes a bound- 
ary layer which features a tangential velocity profile with a positive gradient 
normal to the wall at the wall, a maximum velocity above the wall, before 
adopting a subsequent decrease at greater distances from the wall. In other 
words, the buoyancy imbalance drives flow from left to right adjacent to the 
forcing boundary. In a flow dominated by mechanical forcing (e.g. Case 3), 
the forcing boundary velocity profile is maximum at the wall, and decays 
further from the wall. In this case the wall movement leads and drives the 
flow through the action of viscous shear, carrying the temperature field in 
an almost passive fashion. Hence in this case, where both mechanisms are 
active, the flow will be dominated by natural convection when the boundary 
layer velocities invoked by buoyancy destabilization exceed the mechanical 
forcing boundary velocity. 

A notable feature appearing in this video is an instability in the forcing 
boundary near the right-hand wall that manifests as a short-wavelength 
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(a) Case 1: Ra= 10^°, i?e = 4 072 




{h) Case 2: i?a = 10^°, i?e = 5 700 




Figure 1: Instantaneous snapshots obtained at time 2.4 x 10~^L^/kj' for 
the three cases displayed in the video. The right-hand (heated) half of the 
enclosure is shown, the forcing boundary is at the bottom of each frame, 
and black through red to yellow contours represent increasing temperature. 



4 



waviness that rapidly grows in amplitude and breaks into mushroom-plumes 
rising from the forcing boundary. These are invoked in cases where the 
horizontal convection mechanism has carried cooler fluid over hotter parts 
of the forcing boundary, creating a locally unstable flow with properties 
similar to classical Rayleigh-Benard convection [1], where convection cells 
are established to allow deeper hot fluid to displace cooler fluid at shallower 
depths. 
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